Out-of-phase mixed holographic gratings : a quantative analysis 
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(N ' 

We show, that by performing a simultaneous analysis of the angular dependencies of the ± first 
and the zeroth diffraction orders of mixed holographic gratings, each of the relevant parameters 
can be obtained: the strength of the phase grating and the amplitude grating, respectively, as well 
j /^ ■ as a potential phase between them. Experiments on a pure lithium niobate crystal are used to 

demonstrate the applicability of the analysis. 

<n : 

O . I. INTRODUCTION 

> ' 

Recently, volume holographic gratings with a modulation of both, the absorption coefficient and the refractive 
index, have attracted attention in various materials such as silver halide emulsions [3, 0, Q , doped garnet crystals [3| 
q . or materials with colloidal color centers . The simplest theoretical description (two- wave-coupling theory) of light 
• <-h ' propagation in an isotropic medium with a periodic modulation of the (complex) dielectric constant was already given 
long ago by Kogelnik [||. Considering periodic phase and amplitude modulations, the grating types are treated to be 
in phase. Later Guibelalde generalized the equations to be valid for out-of-phase gratings [7J. The quantity of major 
O ^ interest usually is the (first order) diffraction efficiency rj\ , defined as the ratio of powers between the diffracted beam 
l— ~~ '■ and the incoming beam. For the case of high diffraction efficiencies (above 50%) or even for overmodulated gratings 
[1, H, Q, the so called 'transmission efficiency' 770 , i.e., more correctly termed as zero order diffraction efficiency, was 
also employed for characterization of the grating parameters. It was suggested, that by measuring the diffraction 
^\ ' and transmission efficiency it is possible to evaluate the refractive-index modulation n\ and the absorption constant 
I/"") , modulation a,\ if one assumes in-phase gratings 

■ In this article we show how the shape of the angular dependencies for the ± first and zero order diffraction efficiencies 
\ depend characteristically on the parameters m, ol\ and the phase <p between them. We generalize the formulae given 
Cn ■ in Ref. Q to the case of out-of-phase gratings and demonstrate at two experimental examples that the analysis is 
applicable. This is important, as up to now the evaluation of mixed gratings including a phase was only conducted 
. by beam-coupling experiments, an interferometric technique which is more demanding from an experimental point of 
^— ^ ' view. 
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II. DIFFRACTION EFFICIENCIES OF ZERO AND ± FIRST ORDER 

According to Refs. [1, @| a plane wave propagating in a (thick) medium with a one dimensional periodic complex 
dielectric constant, composed of its real part n{x) = n^ + ni cos (Kx) and imaginary part a(x) = ao + ai cos (Kx + ip), 
yields outgoing complex electric field amplitudes for the (zero order) forward diffracted Ro and (first order) diffracted 
_R ±1 waves. These depend characteristically on the following parameters: the mean absorption constant ao, the 
thickness d of the grating, the dephasing ■d due to the deviation from Bragg's law and the complex coupling constant 
= niir/X — iai/2e ±lv = Ki — in2e ±lv . Further, K denotes the spatial frequency of the grating, no the mean 
refractive index of the medium, and ip a possible phase shift between the refractive-index and absorption grating. 
The goal of an experiment is to extract the grating parameters n\,ai,ip by varying the dephasing, e.g., through 



•Electronic address: martin. fally@univie.ac. at URL: http://nlp.univie.ac.at 



2 



measuring the angular response of 770 = RqRq/I and rj±i = R±iR±i/I where * denotes the complex conjugate and 
I the incident intensity. For simplicity in calculations and as the most often used experimental setup we assume a 
symmetrical geometry, i.e., that the grating vector and the surface normal are mutually perpendicular. A schematic 
of the setup is shown in Figure [T] 




FIG. 1: Schematic of the setup. The angular dependence of the zeroth 770(6) and ± first diffraction orders 77±i(0) from a 
mixed grating is measured by rotating the sample around an axis perpendicular to the grating vector K . Q denotes angles 
outside the medium. Note, that we are within the thick grating regime, where only two beams are propagating simultaneously: 
the zero order together with either the —1 st or the +l st order. 

Slightly adapting the convenient notation from Ref. [l[ the efficiencies for transmission gratings can easily be 
calculated to yield 
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with the abbreviations A(8) = exp{— 2a^D} : D = d/ cosQ and 
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Here, Ob denotes the Bragg angle (inside the medium). Equations (JTJ) and ([2]) are valid for 8 > 0; for 8 < the 
angles and phase-shifts are replaced by their negative values, i.e., »7±i(— 0) = Vti(@) an< ^ ? 7±i( — f) = 7 ?ti(¥')- Note, 
that Equation (fT]) is identical to Equation (11) from Ref. [7j ■ Employing Equations (Q]) and (fj]) we now study the 
particular case of aod = 1 and «2 = ccq/2, i.e., maximal grating strength for the amplitude contribution [||. We 
vary the strength of the phase grating between k\ = K2/4, K2,4«2 with different phase angles tp = 0, 7r/4, 7r/2, 3ir/4 
between the grating types. The angular dependencies of the zero order and ± first order diffraction efficiencies are 
depicted in Figure I(a)-(d). 

At this point let us summarize the main characteristic features occurring in the diffraction efficiencies at the example 
for ip = 7r/4 to obtain a qualitative understanding of the curve shapes and their dependency on the ratio of K\j k^. 

• Zero order diffraction efficiency t\q{8) 



— The curves are symmetric with respect to normal incidence, i.e., 8 = 0. 




e [rad] 
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— Neither the minima nor the maxima of the curve are located at the Bragg angle, except for k\ = or 
k 2 = or tp = it/2. In general the position and the height of the minima or maxima depend in a complex 
way on m, k 2 , tp and even the mean absorption constant cto (see discussion for agd 3> 1). 

— For Ki < k 2 the curve at the Bragg angle extends more to the region above the mean absorption curve (dash- 
dot line, first term in Equations (fT]) and (jJJ ) than below, i.e., as a simple approximation r/™ ax + rj™ ln > 
2A(6b)- The same is true vice versa for K\ > n 2 

— Note, that for \9\ <c 9b the curve resides below the mean absorption curve, for \9\ ^> 9b above 

• Diffraction efficiency 

— The maximum value of the diffraction efficiency differs for t)-i(6b) and t] + i(9b)', in our case »7_i(0b) < 
»7+i(0b)- 

— The curves are symmetric with respect to 9b, i.e., ?7i(#_b + x ) — r h(@B ~ %) except for their different mean 
absorption A(9b ± x). 

— Note, that despite k\ < k 2 the diffraction efficiency 77„i(#b, k 2 /A) > ?/-i(#b,K2) for the minus first 
diffraction order, whereas it is vice versa for the plus first diffraction order, i.e., 7]+i(9b-, K2/4) < ?7+i(0b, k 2 ) 

Next we would like to point out the difference between the curves for various tp values. Figure E{c) shows a unique 
case which is most instructive. For tp = tt/2 the coupling constant n — k\ ± n 2 , G M . Thus a maximum difference 
between ?7_i and rj + i is obtained, culminating in the full depiction of ry_i if k\ = n 2 (sec appendix). Finally, we want 
to draw the attention to the case of <p = 3tt /4 > tt/2. Then r)±i gives identical results as for (p — tt/2. The zero 
order diffraction efficiency 770, however, approaches the mean absorption curve for \9\ 9b from above in the case of 
ip < tt/2 and contrary from below for (p > tt/2 . Considering these arguments it is obvious, that only a simultaneous 
fit of all diffraction data, i.e., zero and ± first order diffracted intensities, allows to extract the decisive parameters 
Ki,K 2 ,tp. On the other hand these curves are therefore fingerprints of the relation between the parameters. The 
following recipe can help in judging about the general situation (for aod ~ 1): 

• Check r)±\: if their magnitudes differ, this is a fingerprint that mixed gratings exist that are out of phase ((/p ^ 0). 
The ratio J]+i/i]_i at the Bragg position obtains a maximum value for tp = tt/2 and for K\ = k 2 [4|. 

• Check ?7o: if 770(6* = 0) < A(0) then \tp\ < tt/2 and else vice versa 

• If r]Q lax + j]™ m > 2A(9b), the absorptive component is dominating and else vice versa. 

• For overmodulated phase gratings another feature of the diffraction efficiencies becomes prominent: the side 
minima near the Bragg peak are lifted to nonzero values (for tp ^ tt/2). This striking feature can already 
be understood in the case tp = where we simply add up the pure absorptive and the pure phase grating. 
The positions of the s th side minima are then given by fii 1 ^ = 2[{stt / D) 2 — n 2 ] 1 ^ 2 (phase grating) and ^i 2 ' = 
2[(stt/D) 2 + a* 2 ;] 1 / 2 (absorption grating). Thus, their minima considerably deviate from each other for n\ t2 w 
stt/D. Recalling, that k 2 = a\/2 < a^/2 such a situation will practically occur if fti 3> tt/D, i.e., for s > 1 
(overmodulated phase gratings exist). This is realized in various systems (see e.g., 0, H, but did not deserve 
proper attention. 

Finally, we would like to recall that for tp — » tp+n the complex coupling parameters fir^ — > k t are interchanged and thus 
the 77-j-i — > 77=pi . For tjq the term in the second line of Equation[2]changes sign because of cos tp — > cos (tp + tt) = — cos tp. 

III. EXPERIMENTAL AND DISCUSSION 

The investigations were performed on a pure congruently melted lithium niobate crystal (thickness: 5mm). Holo- 
graphic transmission gratings were prepared by a standard two-wave mixing setup using an argon-ion laser at a 
recording wavelength of X p = 351 nm. Two plane waves with equal intensities and parallel polarization states 
(s-polarization) were employed as recording beams under a crossing angle of 2Qb = 20.21° (outside the medium) 
corresponding to a grating period of 1000 nm where the polar c-axis is lying in the plane of incidence. The total 
intensity of the writing beams was 9 mW/cm 2 . HPDLC samples were fabricated from a UV curable mixture pre- 
pared from commercially available constituents as previously reported in literature The grating period was 1216 
nm, the grating thickness about 30/xm (llj . After holographic recording we postcured the sample by illuminating it 
homogeneously with one of the UV writing beams. 

The grating characteristics of the samples was analyzed by monitoring the angular dependencies of the ± first and 
zero order diffraction efficiencies. For this purpose the samples were fixed on an accurately controlled rotation stage 
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FIG. 3: Angular dependence of the diffraction efficiencies for the zero and ± first orders of a grating recorded in a pure 
LiNbOa sample. The unequal values of the first order diffraction efficiencies are an impressive signature for the existence of 
mixed phase and amplitude gratings that are out of phase. The zero order diffraction efficiency at the Bragg angles show only 
a slight asymmetry with respect to 9b because the diffraction efficiencies are small (< 1%) and the phase grating is by far 
dominating. The solid lines show a simultaneous fit to r)±i and r/o- The dashed-dot line indicates the mean absorption curve. 



with an accuracy of ±0.001°) and facultatively (HPDLC) in a heating chamber. In the case of LiNbOa we used a 
single considerably reduced readout beam at X r = X p = 351 nm and s-polarization, whereas for the HPDLC a He-Ne 
laser beam at a readout wavelength of A r = 543 nm and p-polarization state was employed. Figure [3] shows the 
experimental curves for the 0.,±1. diffraction orders from a grating recorded in nominally pure congrucntly melted 
LiNbOa . According to the recipe given above we immediately can diagnose mixed out-of-phase refractive-index and 
amplitude gratings, because the > Further by inspecting the zero order diffraction we come to know that the 
phase < (p < it/2. The effects in the zero order are not so prominent for two reasons: the overall diffraction efficiency 
is very small and the phase grating is dominant because the zero order diffraction curve extends mostly to values below 
the mean absorption curve (dash-dot line in Figure [3]). A simultaneous fit of Equations [1] and to the measured data 
yielded the following parameters: n x = (3.01±0.04) x 10~ 6 , a x = 8.18±0.48m -1 , tp = 1.027±0.059, a = 118±1.7m _:l 
with a reduced chi-square value of 1.89 x 10 -7 . From this value and Figurc[3]it is obvious that the equations excellently 
fit the data. 

Finally, we intend to demonstrate the usability of the (qualitative) analysis employing an example with strong 
overmodulation and high extinction: holographic polymer-dispersed liquid crystals (H-PDLCs). Only recently was 
a preliminary beam-coupling analysis of such a system conducted, a task which is not simple from an experimental 
point of view in particular if the experiments should be carried out under high temperatures or application of 
external electric fields. Figure 2] shows the diffraction curves from a grating in a HPDLC at an elevated temperature. 
We can understand the major characteristic features as follows: The liquid crystal (LC) component in an HPDLC 
is highly bircfringent. Statistical alignment of the LC-droplcts of about the light wavelength's size leads to strong 
scattering, i.e., extinction which can be treated similar to absorption provided that multiple scattering does not play 
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FIG. 4: Zeroth and ± first diffraction orders of a strongiy overmodulated grating in HPDLC at 63° Celsius. A = 1.2/itm, A=543 
nm, d = 30 /im (same sample as used for the investigations in Ref. [ly]. The lower graphs show a simulation according 
to Equations [T] and [2] Note, that here the mean extinction is already rather high. We do not expect that a fit could be 
successful for at least three reasons: (1) as in HPDLCs anisotropic gratings are formed, for the basic equations the full theory 
of Montemezzani and Zgonik should be employed [l|. (2) It can be noticed, that around 9 — more than two waves are 
propagating in the medium. Therefore, also the two-wave coupling theory is not fully applicable. Instead a rigorous coupled 
wave analysis should be performed [13). (3) The gratings are expected to be inhomogeneous and non-sinusoidal [O], thus not 
completely fulfilling the requirements 



an essential role. HPDLCs basically consist of alternating regions with high and low concentration of LCs embedded 
in a polymer matrix. Thus, these periodically varying scatterers act as extinction gratings. In addition, of course, 
also the refractive index is strongly modulated (at least via the density changes). Therefore, HPDLCs are typical 
examples of mixed gratings. Furthermore, it is well known in literature that the light-induced refractive-index changes 
are extremely high and strong overmodulation occurs (see e.g. [To!]). Such an example is shown in Figure [U From 
the experimental data we conclude, that combined refractive-index and extinction gratings are produced. This is 
consistent with our previous beam-coupling measurements [l4|. However, we do not dare to decide about a possible 
phase between them. A quantitative evaluation is not possible for this case as we are aware of the fact, that in 
HPDLCs the gratings are anisotropic and thus the basic equations of Ref. [|| should be replaced by the full equations 
given by Montemezzani and Zgonik [l2l ]. In addition, the gratings are usually rather inhomogeneous across the sample 
but might be considerably improved upon further efforts during recording hj. The non-zero minima in the diffracted 
beams partially might originate from overmodulation as discussed above but mainly from the inhomogeneity of the 
gratings and a profile perpendicular to the grating vector [la |. However, a qualitative understanding of the changes 
occuring during heating or applying an electric field can still be read off from the diffraction curves like those shown 
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inRef. [n|. 



IV. REMARKS AND CONCLUSION 



The above discussed analysis is easily applicable for agd r* 1, Ki ks k 2 and 771(6 'b) / b) ^ 0.01, so that with the 
chosen example of LiNbC>3 above we are already at the limit. If one grating type is dominant the analysis still remains 
valid, however, the resulting values for p and the smaller component result in quite large errors. 

We would like to draw the attention to the fact, that for ao <C 1 the absorptive grating strength is considerably 
limited, so that in general the zero order diffraction will not feel the Bragg diffraction. On the other hand, for aod 3> 1, 
the forward diffracted beam will exhibit a maximum near the Bragg position, a fact which is well known in x-ray 
optics (anomalous transmission), see e.g. fl7j . 

We would like to point out, that the analysis of only the first diffraction orders cannot give the full information on 
all relevant parameters 0]. However, it is sufficient to use the ± first together with the zero order diffraction and to 
avoid more demanding beam-coupling (interferomctric) experiments. A prospective phase between the grating and 
the interference pattern(l8l. [l9l. [2fl| . however, cannot be determined by simple diffraction experiments. 

We further would like to emphasize, that the limitations of the coupled wave equations according to Ref. Q should 
be kept in mind when employing Equations [T] and [21 e.g., it is assumed that the gratings are planar, purely sinusoidal 
and isotropic (for anisotropic gratings the theory given in Ref. [l2j should be employed), a\ < ao (for violation of 
this condition see Q) and only two beams are kept in the coupling scheme. If the latter is not applicable the theory of 
rigorous coupled waves has to be applied [l3| , naturally with an increase of the number of coupling constants between 
the beams and thus with loss of simplicity. 
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Appendix 



For the particular case of p — it/ 2 the diffraction efficiencies read: 

r/±i(0) = A(0)z _1 4(«i ±K 2 ) 2 sin 2 [z 1/2 D/2 

= A(9)z- 1 r ± H(K 2 1 - 4) sin 2 (z l / 2 D/2 



z 



A(9)z- 1 \d 2 + A{k\ - 4) cos 2 (V/ 2 £>/2) 



(6) 
(7) 



ti 2 +4(k 2 - k%). 



It's interesting to note, that for this case the diffracted and forward diffracted beams have the functional dependence 
of pure phase gratings with an effective coupling constant of 2[k 2 — k 2 .] 1 / 2 - The amplitude of the diffracted beams, 
however, is enhanced or diminished by a multiplication with or division by r — («i — K2)/{ki + ^2)1 respectively. 
Therefore, it's easy to see that for k\ = K2 the curves shown in Figure [5] (c) arise. 
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